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abstract 

We study ordinary solitons and gap solitons (GSs) in the framework of the one- 
dimensional Gross-Pitaevskii equation (GPE) with a combination of linear and nonlin- 
ear lattice potentials. The main points of the analysis are effects of (in)commensurability 
between the lattices, development of analytical methods, viz., the variational approxi- 
mation (VA) for narrow ordinary solitons, and various forms of the averaging method 
for broad solitons of both types, and also the study of mobility of the solitons. Under 
the direct commensurability (equal periods of the lattices, Lun = Lnoniin), the family of 
ordinary solitons is similar to its counterpart in the GPE without external potentials. 
In the case of the subharmonic commensurability, with Lun = (l/2)Lnoniin, or incom- 
mensurability, there is an existence threshold for the ordinary solitons, and the scaling 
relation between their amplitude and width is different from that in the absence of the 
potentials. GS families demonstrate a bistability, unless the direct commensurability 
takes place. Specific scaling relations are found for them too. Ordinary solitons can 
be readily set in motion by kicking. GSs are mobile too, featuring inelastic collisions. 
The analytical approximations are shown to be quite accurate, predicting correct scal- 
ing relations for the soliton families in different cases. The stability of the ordinary 
solitons is fully determined by the VK (Vakhitov-Kolokolov) criterion, i.e., a negative 
slope in the dependence between the solitons's chemical potential ^ and norm A'^. The 
stability of GS families obeys an inverted ("anti-VK") criterion, d/i/dN > 0, which is 
explained by the approximation based on the averaging method. The present system 
provides for a unique possibility to check the anti-VK criterion, as fi{N) dependences 
for GSs feature turning points, except for the case of the direct commensurability. 
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1 Introduction 



It is well known that periodic potentials, induced by optical lattices (OLs), provide 
for a versatile tool for controlling dynamics of Bose- Einstein condensates (BECs). 
This tool is especially efficient for the creation and stabilization of solitons - both 
ordinary ones and gap solitons (GSs), which are supported by the interplay of the 
OL potential and self-repulsive nonlinearity. Many results obtained in theoretical 
and experimental studies of this topic were summarized in reviews focused on one- 
dimensional (ID) [T] and multidimensional [2] matter-wave dynamics. Earlier, a similar 
model was introduced in optics for the description of spatial solitons in nonlinear 
waveguides with a periodic transverse modulation of the refractive index [^. In the 
experiment, lattices controlling the transmission of optical beams were implemented in 
the form of photoinduced gratings in photorefractive crystals, which made it possible 
to create various species of ID and 2D spatial solitons [4]. Gratings were also created 
as permanent structures written by femtosecond laser beams in silica Recently, 
this topic was reviewed in Ref. [B]; a closely related topic is the study of discrete 
solitons in optics, which was a subject of another comprehensive review in Ref. [7|- 

A different possibility, which has drawn much attention in studies of BEC (thus far, 
primarily at the theoretical level) , is the use of a spatially profiled effective nonlinearity, 
that may be implemented by means of properly designed configurations of external 
fields via the Feshbach-resonance effect. In terms of the condensed-matter theory, 
the nonuniform nonlinearity coefficient induces an effective pseudopotential [SlIS], that 
can be used to control the dynamics of localized modes. Various problems of this sort 
were considered in ID settings, with periodic pseudopotentials in the form of nonlinear 
lattices (NLs) [101 llll I12| . as well as with spatial modulations of the nonlinearity 
coefficient represented by one [13] or two [9] delta-functions (actually, the model with 
the self-attractive nonlinearity concentrated at a single delta-function was introduced 
long ago as a model for tunneling of interacting particles through a junction 1141). 
In particular, the configuration with two delta-functions makes it possible to study 
spontaneous symmetry breaking of matter waves trapped by a symmetric double-well 
pseudopotential [9,. Specially chosen profiles of the nonlinearity coefficient may also 
be employed to design a pulse-generating atomic- wave laser [T^, as well as traps and 
barriers for such pulses ^16j. The analysis of ID matter- wave solitons in NLs was 
further extended for two-component models |17l I18| . and for some spatiotemporal 
patterns of the nonlinearity modulation [19]. Certain results for solitons supported 
by 2D nonlinear pseudopotentials were reported too, although the stabilization of 2D 
solitons in this setting is a tricky problem PUI . 

In addition to the BEC, a periodic modulation of the nonlinearity is possible in 
optics, where it was analyzed in terms of temporal [21] and spatial [221123] solitons. A 
discussion of practical possibilities to create NLs in optical media in the "pure form" 
(without affecting the linear properties, i.e., the refractive index) can be found in 
recent paper [18]. An experimental observation of NL-supported optical solitons (in 
the form of surface solitons at an interface between lattices) was reported in Ref. . 

A natural generalization of the study of NL pseudopotentials is to consider atomic 
and optical media equipped with a combination of nonlinear and linear lattices, in 
ID [251 1261 112] and 2D [27] settings. In the experiment, this may be implemented 
by applying the above-mentioned techniques simultaneously - for instance, combining 
the OL, which induces the linear periodic potential in the BEC, and the patterned 
magnetic field, which gives rise to the nonlinear pseudopotential via the Feshbach 
resonance. Actually, photonic-crystal fibers, where various species of spatial solitons 
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have been predicted [28], also belong to this type of media. Under special conditions, 
the combined models admit analytical solutions, see Ref. [12] and references therein. 
In particular, exact solutions were elaborated in detail for the lattices of the Kronig- 
Penney (piecewise-constant) type [23]. 

The objective of this work is to investigate the existence, stability, and mobility of 
ordinary ID solitons and their GS counterparts in the combined NL-OL model, with 
emphasis on effects of the commensurahility (spatial resonance) between the nonlinear 
and linear lattices. In the general case, the solitons are found in a numerical form. For 
narrow ordinary solitons (those whose chemical potential, /i, falls into the semi-infinite 
gap of the linearized version of the model) , we also develop a variational approximation 
(VA). For broad solitons, both ordinary ones and GSs with ^ belonging to the first 
finite bandgap, we elaborate analytical approximations based on different versions of 
the averaging technique. 

The stability of the solitons is investigated below by means of systematic simula- 
tions of their perturbed evolution. We conclude that the well-known VK (Vakhitov- 
Kolokolov) criterion 29; completely determines the actual stability of the ordinary soli- 
tons (the criterion states that a necessary stability condition is a negative slope in the 
dependence of the solitons' chemical potential, ^, on their norm. A*', i.e., d^/dN < 0). 
For the GS families, our analysis leads to a different conclusion: their stability fully 
obeys an "anti-VK" criterion, viz., dfi/dN > 0. As a matter of fact, all stable GS fam- 
ilies in previously studied models had only the positive slope of ij.{N), thus satisfying 
the latter criterion automatically. However, the present model offers a rather unique 
chance to test it in a nontrivial situation, when ^j.{N) curves for the GSs may have 
portions with both positive and negative slope, separated by turning points. Using 
the averaging method, we also produce a justification for the anti-VK criterion, which 
is relevant, at least, for subfamilies of broad GSs. 

The rest of the paper is organized as follows. The model is formulated in the next 
section, which also reports basic numerical results. In the case of the direct com- 
mensurahility between the nonlinear and linear lattices (for those with equal periods, 
iiin = i'noniin), it is coucluded that ordinary solitons are similar to soliton solutions 
of the nonlinear Schrodinger equation (NLSE) in the free space (without an external 
potential), in the sense that there is no threshold for their existence, the entire soliton 
family is stable, and the amplitude and width of the soliton, A and W, obey the usual 
scaling relation, A oc 1/14^. On the contrary to that, there is an existence threshold 
for ordinary solitons in the case of the incommensurability, or if the commensurability 
(spatial resonance) between the lattices is subharmonic, with Lu^ = inoniin/2. In those 
cases, the scaling relation between the amplitude and width is different too - featuring, 
in particular, A oc for the subharmonic resonance. The same commensurability- 

depending change in the scaling is observed in GS families, which also demonstrate a 
bistability in cases different from the direct commensurability. 

Analytical results, which may explain a considerable part of the numerical findings, 
are collected in Section III. First, we develop the VA for narrow ordinary solitons, and 
then the averaging method is reported, in different forms for different cases. It is 
demonstrated that both the VA and averaging produce results which are in good 
agreement with numerical observations, in relevant regions of the parameter space. In 
particular, as briefiy mentioned above, the averaging lends an explanation to the "anti- 
VK" stability criterion for GSs. In the case of Lun = I/noniin/2, an average equation 
for the envelope amplitude includes a quintic nonlinear term, rather than the cubic 
one, which accounts for the above-mentioned change in the scaling relation between 
A and W. The paper is concluded by Section IV. 
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2 Numerical results 



2.1 The model 

The model combining the hnear OL potential and nonhnear NL pseudopotential is 
based on the known effectively one-dimensional Gross-Pitaevskii equation for the 
mean-field wave function, (f){x,t) ,23, 25, 26, 12 . In the scaled form, the equation 
is 

i(j)t = - (1/2) - [ecos(27ra::) + gcos{-Kqx) \(t)\^] (j), (1) 
where e is the strength of the linear OL, and g — ±1 is fixed by the normalization. 
The NL wavenumber is g, the above-mentioned periods of the linear and nonlinear 
lattices being Lun = 1 and Lnoniin = 2/q. For q = and g = —1, Eq. ([!]) amounts to 
the NLSE with the OL potential and constant coefficient in front of the self-defocusing 
cubic term. As is well known, this equation supports families of stable GS solutions 
[Tl|30]. On the other hand, in the absence of the OL, e = 0, the NL supports ordinary 
solitons, but not GSs [inilllllll]. 

Stationary solutions with chemical potential fi are looked for in the form of (j}{x, t) = 
u{x) exp(— i/xJ), where u{x) is a real function. The substitution of this expression into 
Eq. ([T]) yields an ordinary differential equation, 

fiu = ~ (1/2) u" — [ecos(27ra;) 4- g cos(7rga;)u^] it, (2) 

which can be derived from the corresponding Lagrangian, 

Uu^-(l/2)(u') -f-ecos(27ra;) -f (g/2) cos (vrga;) u\dx. (3) 

-oo 

We have constructed numerical solutions for localized stationary modes by means of 
the shooting method applied to Eq. The stability of the so found solutions against 
small perturbations was tested through direct simulations of Eq. ((!]). The results were 
also compared to predictions of the VK criterion for the ordinary solitons, and to the 
above-mentioned "anti-VK" criterion for GSs. Numerical results reported below are 
obtained for the OL strength e = 5, which adequately represents the generic situation, 
for the ordinary solitons and GSs alike. 

Figures [TJa) and (b) display, respectively, examples of ordinary solitons (for g = 
-1-1, q = 1) and GSs (for g — —1, q — 1) with equal amplitudes. The choice of q — 1 
corresponds to the above-mentioned case of the subharmonic resonance between the 
OL and NL, iun = (l/2)I/noniin- In this case, stable ordinary solitons and GSs coexist 
for either sign of g. If, for instance, g = -1-1, ordinary solitons are located around even 
sites of the NL, x = 2n (with integer n), while GSs may be centered at odd sites, 
X = 2n + l. 

2.2 Ordinary solitons 

Figure ^sl) represents a family of the ordinary solitons, by means of the relation be- 
tween their norm, TV = v?{x)dx, and chemical potential /i, at three characteristic 

values of the NL's wavenumber, q = 1, \/5 — 1, and 2, for g — +1. These values are cho- 
sen because q = 2 corresponds to the direct commensurability between the linear and 
nonlinear lattices (Lun = inoniin), 9=1 represents, as said above, the subharmonic 
commensurability [Liin = (1/2) Lnoniin], and q = y/5 — 1 corresponds to incommensu- 
rate lattices. All values of fj, for the ordinary solitons fall into the semi-infinite gap of 
the spectrum induced by the OL potential in the linearized version of Eq. ((l|. 
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Figure 1: Typical examples of a stable broad ordinary soliton (a) and stable 
gap soliton (b), at g = +1 and —1, respectively, in the case of the nonlinear 
lattice with period inoniin = 2 (i.e., q = 1). The amplitudes of both solitons are 
A = 1.5. The chemical potential and norm are, respectively, fi = —1.198 and 
N — 1.628 for the ordinary soliton, and fi = 2.629 and = 2.176 for the gap 
soliton. 



At g = 2, direct simulations demonstrate that the entire soliton family is stable, 
precisely as suggested by the VK criterion, d^/dN < 0, see Fig. [IJa). For q = 1 
and q — y/5 — 1, the results are different, featuring non-monotonous relations fJ,{N). 
Accordingly (again in the agreement with the prediction of the VK criterion), nar- 
row solitons with larger amplitudes, which correspond to the branches of jJ,{N) with 
dfi/dN < in Fig. [SJ^a), are stable, while their loosely bound (broad) counterparts, 
corresponding to the branches with dfi/dN > 0, are unstable. Another difference from 
the case of the direct commensurability is that there is a minimum value of the norm, 
A^min (the threshold), which is necessary for the existence of the ordinary solitons at 
g = 1 and — 1, while there is no threshold at g = 2. 

In fact, the situation in the case of g = 1 and — 1 - the existence of the 
threshold value, A'min, which separates stable and unstable branches of the ordinary- 
soUton solutions - is qualitatively similar to what is known in the model with the NL 
but no linear potential ^lOj. On the other hand, the situation in the case of g = 2 
- the absence of A'min and the existence of the single branch of the soliton solutions, 
which is entirely stable - resembles well-known properties of soliton solutions of the 
NLSE without any lattice, linear or nonlinear. 

The evolution of those ordinary solitons which are unstable is illustrated in Fig. 
[3] for g = 1 and g = -1-1. It is observed that the unstable soliton, with initial ampli- 
tude A — 0.9, rearranges itself into a narrower persistent breather, with time-average 
amplitude Ahr ~ 1.1, while the norm is kept constant. The transformation of the 
unstable ordinary solitons into stable breathers is also similar to what was reported in 
the model without the OL [10) . 

Properties of the ordinary solitons in the same three families, with g = 1, — 1 
and 2 (and g — +1), are further illustrated in Fig.[5Jb) through relations between the 
soliton's amplitude. A, and its width, W, which we define by 

/ + 00 
\<t){x)f{x-L/2fdx (4) 
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(a) (b) 




Figure 2: (a) Chemical potential /i versus norm N for families of ordinary soli- 
tons at different values of q and g = +1. Portions of the curves with d/.i/dN < 
and dii/dN > represent, respectively, stable and unstable (sub)families, in 
agreement with the VK criterion, (b) The log-log plot of amplitude A versus 
width W [the latter is defined as per Eq. (|3])]. Bold dashed lines correspond to 
unstable branches of the soliton families, ior q — 1 and V^— 1. Thin dashed ref- 
erence lines designate scalings which different families obey at small values of A, 
namely, W (x 1/A and W cx 1/A'^. (c) The stability boundary for the ordinary 
solitons, which is defined, as per the VK criterion, by condition dN/dji = 0. 
As predicted by the criterion and verified in direct simulations, the solitons are 
stable above the boundary. 




Figure 3: An example of the spontaneous rearrangement of an unstable ordinary 
soliton, with amplitude A — 0.9 and norm N = 1.32, into a robust breather, 
at q = 1 and g — +1. (a) The evolution of the soliton's amplitude, (b) The 
field profile, \(j){x,t)l at t = 110. (c) The profile at t = 140. Panels (b) and (c) 
display the shape of the breather at points where its width is close, respectively, 
to the minimum and maximum values. 
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(a; = I//2 is the central point of the integration domain). At g = 2, relation W^(^) 
features scaling W oc 1/A for relatively small values of A. This is the same scaling as 
featured by exact soliton solutions of the NLSE in the free space, which is in line with 
the above observation that the soliton family at 5 = 2 is similar to that in the NLSE 
without any lattice. However, the scaling is different in the other families, featuring 
W{A) oc l/A^ for g = 1, and W{A) oc * for g = \/5 - 1. The two latter scaling 
relations imply that the width of the respective solitons is essentially larger than in 
their free-space counterparts, therefore we call them broad solitons. 

Figure OJc) summarizes the results by means of the stability boundary for the 
ordinary solitons in the plane of {q,N). The boundary is identified as a VK-critical 
curve, along which dfi/dN vanishes, the stability area (with dfi/dN < 0) being located 
above the curve. Systematic simulations performed in regions below and above the 
boundary have confirmed that the solitons in these regions are, respectively, unstable 
and stable (unstable solitons transform themselves into breathers, as shown in Fig. |3]). 
The stability area reaches the limit of A'' = (very broad solitons with a vanishingly 
small amplitude) in the form of the cusps in Fig. [2fc) at g = and q = 2. Recall 
that q = with g = +1 corresponds to the constant coefficient of the self-attractive 
nonlinearity in Eq. ([T]), while q = 2 corresponds to the direct commensurability 
between the linear and nonlinear lattices. 

2.3 Static gap solitons 

In this work, the consideration of GS families was confined to the first finite bandgap 
induced by the OL potential, in terms of the linearized version of Eq. lU. The N{^) 
and Vy(A) curves for these families are displayed in Figs. HI a) and (b), for the same 
three cases as above, viz., q — 1, \/5 — 1 and 2, fixing g — —1 [width W is again defined 
as per Eq. ([7|]. The VK criterion does not apply to GSs. Nevertheless, the results 
strongly suggest that the stability of all GS families follows an "anti-VK" condition, 
dn/dN > 0. For subfamilies of broad GSs, this condition is derived below from an 
effective envelope equation for broad GSs which amounts to an "inverted" NLSE, with 
the self-repulsive nonlinear term and negative effective mass, see Eqs. (|15|l . In that 
approximation, GSs reduce to ordinary solitons if the wave function is subjected to 
the complex conjugation, which implies the inversion of the sign of ^, and of dji/dN 
as well. 

In accordance with what is said above, the entire GS family for q = 2, which 
satisfies the "anti-VK" condition everywhere in Fig. |4ja), is found to be completely 
stable in direct simulations. On the other hand, the fJ,{N) curves for q — 1 and 
q = \/5 — 1 feature two folds, and direct simulations corroborate the instability of 
portions of the GS families with dfi/dN < 0. To the best of our knowledge, the 
present model produces the first example of I-l{N) characteristics for GSs with turning 
points, which makes the anti-VK criterion amenable to the actual verification. In 
the standard model with the constant nonlinearity coefficient [321 [l], as well as in 
its version with the quasiperiodic OL potential i33j, the monotonous character of the 
curves does not allow the verification of the criterion. 

The stability diagram in the plane of {q,N), as predicted by the "anti-VK" cri- 
terion, is displayed in Fig. HJc). It includes two critical curves with dN/d/i = 0. As 
suggested by Fig. IHb), and completely confirmed by systematic simulations, in the 
regions above the top curve and below the bottom one there is a single solution, which 
is stable. Between the curves, there are three solutions, two of which are stable, i.e., 
this is a bistability region. The critical curves feature cusps near q = G and q = 2, 
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(a) (b) 




Figure 4: (a) The relation between norm N and chemical potential fi of gap 
solitons for g = 1,-^5—1, and 2, for g = —1. Portions of the curves with 
dji/dN > and dji/dN < represent stable and unstable (sub)families, re- 
spectively, obeying the "anti-VK" criterion (see the text), (b) The log- log plot 
of amplitude A versus width W, for the gap solitons. Bold dashed portions of 
the curves correspond to unstable solutions, with dfx/dN < 0. Thin dashed 
reference lines designate scalings W (x 1/A and W oc 1/A'^. (c) Critical lines 
dN/dfi = in the plane of {q, N). There is a single stable gap soliton above the 
upper line and beneath the lower one, and three solutions - two stable and one 
unstable - in the bistability region between the two lines. 



similar to the situation displayed in Fig. |2jc). 

From Fig. 4(b) we conclude that the scaling relations between the GS's width and 
amplitude for broad solitons (with small values of A) take the following form: For 
g = 2, W ~ l/A; for g = %/5 - 1, W oc and for q = 1, oc i.e., 

almost exactly the same as their counterparts for the ordinary solitons, see above. All 
portions of the GS families obeying these scaling relations are stable. On the other 
hand, the simulations demonstrate that unstable GSs (those with dfi/dN < 0) are 
not transformed into breathers, unlike unstable ordinary solitons, but rather suffer a 
gradual decay into quasi-linear waves (not shown here). 



3 Analytical methods 

3.1 The variational approximation for ordinary solitons 

Narrow stationary solitons of the ordinary type (corresponding to g — +1) can be 
naturally approximated by means of the VA, using the simplest Gaussian ansatz [3], 

u{x) ^ Aexp[-x^/ {2W^)] , (5) 

with norm A'^ = ^Jt^A^W. The substitution of ansatz ([5]) in Lagrangian Q yields the 
effective Lagrangian, written in terms of the norm, instead of amplitude A: 

or Ar 1 AT -■n-^M'^ , A'^ -Tr^g^W^/S , a\ 

2Lcfi = ^iN ~ — — + eNe H — — e ^ ' . (6) 

4M/^ 2y^W 
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Figure 5: Comparison of the variational (dashed hned) and numericaUy found 
(chains of symbols) curves n{N) for ordinary solitons {g — +1). 

Variational equations following from Eq. dLcd/dW — and dLcd/dN = 0, take 
the form of 

4^2eVKV-'^' + (l/^/2^) (l + ^^fc^i^^) We-(^'^)'/« = 1, (7) 

(W^)"'-ee-<"'^''- (l/%/2^)(iV/Ty)e-(^'^)'/« = M- (8) 

Figure 5 compares the ^(A'') curves produced by Eqs. ([7|l and ([HI) to their numeri- 
cally found counterparts, for the NL wavenumbers q — 1 and 2. It is seen that the VA 
provides for a good approximation for narrow solitons with values of /i which are not 
too close to the edge of the semi-infinite gap. Portions of the curves corresponding to 
broad solitons, which are located near the gap's edge, are not captured by the VA, as 
the actual shape of these solitons is different from the Gaussian, see, e.g., Fig. [TJa). 



3.2 The averaging method 
3.2.1 Ordinary solitons 

The approximation based on averaging can be applied to broad solitons, which have 
a small amplitude and large norm. In the case of the ordinary solitons (g — +1), 
this is the situation opposite to that (narrow localized modes) for which the VA was 
presented in the previous section. To develop the averaging approach for ordinary 
solitons, we adopt the ansatz 

(l){x,t) = <i>{x,t)[l + 2a cos(2ttx)], (9) 

where the slowly varying amplitude function, multiplies the simplest approximation 
for the Bloch wave function which may be used near the edge of the semi-infinite gap, 
with Q — ^x/tt^ -f 6^/2 ~ ""^^ (this approximation is obtained by dint of the analysis 
presented in Ref. [lU]). The substitution of ansatz © into Eq. ^ and averaging, 
also performed along the lines of Ref. [30] , lead to the asymptotic NLSE for the slowly 
varying envelope function, 

= (10) 
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(a) (b) 
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Figure 6: (a) The numerical obtained profile of an ordinary soliton with am- 
plitude ^ = 0.1 for g = +1 and q = 2. (b) The profile predicted for the same 
soliton by means of the averaging method, see Eqs. ©, P^ . (fTT|) . and P^ . 



where the calculations yield the following coefficients: 



m 



(orf) _ 2n^ + e^+ 7rV47r4 + 2e^ ^^^^ 



off 



IOtt^ + £2 _ Stt V47r4 + 2e2 ' 



(ord) _ ([1 + 2q cos(27rx)]'^cos(g7ra:)) 
ffcff - ^^2a2 ' ^^^^ 

with (...) standing for the spatial average. Obviously, the effective nonlinearity coeffi- 
cient given by expression (|12p vanishes unless q takes values 0, 2, 4, 6, or 8. In particu- 
lar, 5^°"^' = -(l + l2a^ + 6a^)/{l + 2a^) for q = 0, and = -4Q(l + 3a2)/(l + 2a2) 
for q = 2. Actually, p^g'*^ does not vanish at g = 2n with any integer n, if higher- 
order harmonics are kept in the expansion of the Bloch function at the edge of the 
semi-infinite gap, 1 + '2'}2„ cos{2-nnx) , cf. the lowest-order approximation used in 
Eq. @. With the value of the OL strength adopted in the numerical simulations 
reported above, e = 5, Eq. (|lip yields m^^'^^ « 1.128, which is virtually identical to 
the numerically found effective mass, see Eq. (??). 

The soliton solution to Eq. (I10|) with an arbitrary amplitude. A, is 



$ = Aexp ^^gir'Ahj sech i^yj gi^' ^iT' Ax j . (13) 

This solution explains scaling W oc I /A, which is observed in Fig. [^b) for broad 
ordinary solitons in the case of g = 2. Figure 6 displays a direct comparison of profiles 
of a typical broad ordinary soliton, as obtained in the numerical form and produced 
by the averaging method. Good agreement between the two profiles is obvious [note 
that the figure displays the full wave function, |<?I>(a;)|, rather than the envelope, "l>(a;), 
seeEq. ©]. 
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3.3 Gap solitons 



3.3.1 Direct commensurability between the nonlinear and linear lat- 
tices 

To apply the averaging approximation to GSs, which correspond to g = — 1 in Eq. ([1]), 
we follow Ref. [30] and adopt the simplest ansatz which is relevant in this case, 

(p {x,t) = $ {x,t) cos{nx), (14) 

with a slowly varying amplitude "I> {x, t). The difference in the form of the carrier wave 
in this expression, in comparison to Eq. (jOj, is due to the fact that, near the edge 
of the first finite bandgap, the Bloch function is close to a periodic function, whose 
period is twice as large as that of the underlying OL potential. On the contrary, near 
the edge of the semi-infinite gap the period of the Bloch function coincides with the 
OL's period. 

The substitution of ansatz (|14[) in Eq. ^ and the application of the averaging 
method yields the respective asymptotic NLSE equation for the amplitude function, 

(gap) 2 

dt - 2m(^^P' dx^ ^^^^ 

cff 

cf. Eq. I|10[l. The effective mass and interaction coefficients in Eq. (|15[) are found to 
be 

m(r'=e/(.-2.^), (16) 

<7^|"P> = 2 (cos(7rx))^ cos(g7ra:)) . (17) 

Note that coefficient p7p is difi^erent from zero only for three values of q, viz., g^.^^^ {q = 
0) = 3/4, gi%''''\q = 2) = 1/2, and gil'"-'\q = 4) = 1/8. Together with 



,(gap) 



-0.339, (18) 



which Eq. (|16|l yields for e — 5 (recall this value of the OL strength is fixed in the 
present work) , the complex conjugate form of Eq. I|15p gives rise to the usual soliton 
solutions for $*, cf. solutions H13p for <I>. This fact gives the explanation to the scaling 
W oc 1/A for the broad GSs, as observed in Fig. IH^b) for q — 2. For stationary 
solutions, the complex conjugation implies, as mentioned above, the reversal of the 
sign of the chemical potential. This explains why the stability of the broad GSs obeys 
the "anti-VK" criterion, d^/dN > 0, which is simply the reverse of the ordinary 
negative-slope VK condition for the stable solutions of the equation for $*. 

The above description is also relevant for moving GSs. In particular, the compar- 
ison of the analytically predicted effective mass ()18p to the empirical dynamical mass 
(??), drawn for the moving soliton from numerical data at g = 2, clearly demonstrates 
the high accuracy of the averaging approximation for the broad GSs, both static and 
moving ones. 



3.3.2 Subharmonic commensurability between the nonlinear and lin- 
ear lattices 

The above approximation for GSs does not produce any definite result for q = 1, when 
the period of the NL in Eq. ((TJ is twice as large as OL's period (the subharmonic 
resonance between the nonlinear and linear lattices, as defined above). To derive an 
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effective envelope equation in this case, we notice that the substitution of original 
ansatz (|14[l in Eq. ^ and making use of the same effective mass as given by Eq. (|16|l. 
but without averaging the nonlinear term, gives rise to the following equation: 

cff 

where it is taken into regard that q = 1. This equation suggests that ansatz p4l) 
should be replaced by the following one, for stationary solutions: 

{x, t) = e"'''* [<E>i (x) cos(7rs) + <1>4 {x) cos" (nx)] , (20) 

where both functions $i and $4 are slowly varying ones. The substitution of ansatz 
(|2U|) into Eq. ((T]) and straightforward trigonometric expansions make it possible to 
eliminate $4 in favor of <l>i: 



<f4(x)= (m--^^ I mx^-fiix). (21) 

^'"cff 

The remaining equation for "I'i(3;) is the NLSE with the quintic self-focusing nonlinear 
term: 

M'&i = -, .^r^ 7 — ^ $1 ^1- (22) 

An obvious soliton solution to Eq. (|22|) with arbitrary amplitude A is 

$1 =yl Vsech(A;x), fc = ^/5m^f (tt" - 2m(|"PV) /i = - (smj^f p') fc^ 

(23) 

Solution (|23[) yields the scaling relation W oc 1/^^, which explains the same scaling 
that was observed, at g = 1, for broad GSs in Fig. IHb). 

It is necessary to mention that the nonstationary version of Eq. (|22[) . with re- 
placed by id^i/dt, corresponds to the one-dimensional NLSE with the critical (quintic, 
in the ID case) self-focusing nonlinearity, whose solutions, tantamount to those given 
by Eq. (|23p , are the so-called "one-dimensional Townes solitons" . It is well known that 
the entire family of such solitons is unstable (see, e.g., Ref. [34]). Nevertheless, simu- 
lations of Eq. lU demonstrate that the GSs approximated by the asymptotic solution 
(|23|) form a stable family, as long as the solitons remain broad, see the corresponding 
stable branch in Fig. |4jb). Thus, the asymptotic description of the broad GSs by 
means of Eq. H22p is valid, at g = 1, only for static solutions, while their dynamical 
behavior does not obey the straightforward time-dependent version of this equation. 



4 Conclusion 

We have investigated the existence, stability and mobility of ordinary solitons and 
GSs (gap solitons) in the ID model combining nonlinear and linear periodic lattices in 
the Gross-Pitaevskii equation for the nearly one-dimensional BEG. The emphasis was 
made on the study of effects of the commensurability and incommensurability between 
the lattices, as well as on the development of analytical methods - the VA (variational 
approximation) for narrow ordinary solitons, and averaging method for broad solitons 
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of both types. We have demonstrated that, in the case of the direct commensura- 
bility between the lattices (equal periods), the ordinary solitons are similar to their 
counterparts in the free space. In the case of the subharmonic commensurability and 
incommensurability, the situation is different, featuring the existence threshold for the 
solitons, and a different scaling relation between their amplitude and width. Simi- 
lar scaling relations are found for GS families, which demonstrate the bistability in 
cases different from the direct commensurability. Ordinary solitons may travel long 
distances, if kicked. Broad GSs are mobile too, although collisions between them are 
inelastic. 

The analytical approximations demonstrate good accuracy in appropriate param- 
eter regions. In particular, they correctly explain different scaling relations for the 
soliton families at different commensurability orders. As concerns the stability of the 
ordinary solitons, it is accurately predicted by the VK criterion. Simultaneously, the 
stability of GSs obeys the "anti-VK" criterion, an explanation to which was given by 
means of the effective equation produced through the averaging method. A notable 
feature of the present model is that it gives rise to characteristics fJ.{N) for the GSs 
that feature turning points (except for the case of the direct commensurability). Un- 
like previously studied models with the constant nonlinearity coefficient, the presence 
of the turning points has made it possible to test the anti-VK stability criterion for 
the GS families. 

The analysis reported in this work may be naturally extended by studying GSs 
in higher finite bandgaps; in particular, a challenging issue is to verify the "anti- 
VK" stability criterion in the higher gaps. It may also be interesting to develop a 

systematic analysis of the commensurability for solitons and solitary vortices in the 
2D model based on the combination of linear and nonlinear lattices. 
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